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(2J[)' The time evolution of wave packets in a harmonic oscillator potential is 
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<^ ■ studied. Some new results for the most general case are obtained. A natural 

number, called "degree of rigidity" , is introduced to describe qualitatively how 
much the shape of a wave packet is changed with time. Two classes of wave 
packets with an arbitrarily given degree of rigidity are presented. 



> 

o 



I. INTRODUCTION 



OO ' PACS numbers: 03.65.Ta 

O ' 

<N 
O 

S3 
Pk 
^— > 
c . 

In recent years, the time evolution of localized wave packets of various quantum-mechanical 
q-( systems has been widely discussed in the literature (see Ref. |] for a review). As for wave 
packets of the harmonic oscillator, the subject is as old as quantum mechanics itself. Though 
some aspects of the Gaussian wave packets (displaced or squeezed ground state) are still being 
discussed in the current literature 0, [|, the problem was essentially solved in the very early 
years ||, j|, |J. In the 1950's, several authors studied more general displaced number states and 
found that they keep their shape unchanged while their center oscillates like a classical particle 
|7|, [|, 0, |IOfl . Displaced and squeezed number states of special forms were also considered in 
some of these papers || [!|; it turns out that their width is also oscillating, so the shapes of 
such wave packets change with time apparently. Displaced and squeezed number states of the 
most general form were studied in a recent paper ]TT|. This may represent the most general 



case where the time dependence of the wave packets could be worked out explicitly. 

The questions we are concerned with are: First, if the time dependence of the wave packet 
cannot be obtained explicitly, what conclusions could be made in regard to its time evolution? 
Second, if the shape of the wave packet changes with time, how can we describe the level of 
changes for different cases? For the first question, it is known that the center of the wave 
packet moves like a classical particle while its width pulsates. However, one can find wave 
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packets whose width keeps unchanged but whose shape still changes with time. To describe 
the change of the shape in the general case we introduce general-order moments of x about 
its center, defined as the mean value of (x — x t ) K and denoted by Qx{t), where x% is the 
mean value of x (which represents the center of the wave packet) and K a natural number. 
Differential equations for these moments are established for general K, and are solved for the 
most general initial conditions up to K = 4. They are essentially oscillating, but higher-order 
ones involve higher frequencies and thus their time dependence is more complicated. Wave 
packets could be found whose Qk{^) are time independent up to K = 4. Obviously the shape 
of such wave packets changes with time less than those that only keep their width [represented 
by a/ QAfj] unchanged. At this stage the answer to the second question is becoming clear. 
We introduce a natural number called the "degree of rigidity" for a wave packet: If Qxit) 
is time independent for K — 2, 3, . . . , 2N [note that Qo(t) = 1 and Q\(t) = 0] but not for 
K = 2N + 2 (for K = 2N + 1 it may or may not be time independent), we say that the 
wave packet has degree of rigidity N. It is obvious that the larger N is, the less the shape 
changes with time. The shape of a wave packet with a time- dependent width changes with 
time apparently and thus has no rigidity. One with a constant width but with all higher-order 
moments time dependent has the "ground" degree of rigidity 1. On the other hand, if Qk{^) 
is time independent for all K, then the degree of rigidity for such a wave packet is infinity. In 
other words, its shape is perfectly rigid. Typical examples of such wave packets are displaced 
number states. Two classes of wave packets with an arbitrarily given degree of rigidity are 
presented. 



II. SOME GENERAL RESULTS 

Consider localized wave packets of the harmonic oscillator in one dimension, whose time 
evolution is governed by the Hamiltonian 

The normalized initial wave function ip(x, 0) at t — will be denoted by i/jq for convenience. 
Localization means that the mean value of x k p l in ip is finite for any nonnegative k and /. 
The wave function ip(x, t) at time t will be denoted by ip t - The mean values of x and p at time 
t are 

x t = (^t,x^ t ), Pt = (tpuP^t)- (2) 
Using the Schrodinger equation it is easy to show that 

x t =pt/^, p t = -fxuj 2 xt. (3) 

The solution of these equations is 

x t = xq cos out + — sin out, p t = p cos out — uujxq sin cut. (4) 
fiou 

The solution means that the center of an arbitrary wave packet moves like a classical particle. 
This is a well known result. 

To study the change with time of the shape of the wave packet we consider the mean value 
of (x — x t ) K in ip t where K is a natural number. These mean values are denoted by Qxif) in 
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the following. If the Qk^Ys are time independent for all K, then the wave packet obviously 
keeps its shape unchanged while it is moving, and this is the case for displaced number states 
|7|, §, |[| [llj. In order to examine their time dependence for a general wave packet, we define 
the real quantities 

Rki(t) = ~(i> t ,{(x-x t ) k ,(p-p t ) l }i; t ), S^t) = ^ t ,[(x-x t ) k ,{p-p t ) l ]^ t ), (5) 

where k and I are nonnegative integers, and {F, G} = FG + GF, [F, G] = FG — GF for any 
operator F and G. Obviously i?jco(t) = Qi<(t). In the following we also denote Roxit) by 
Pk(€). It is straightforward to show that they satisfy the following equations 

R kl (t) = - lfiu 2 R k+1 ^ 1 (t) + —k(k - l)S k - 2jl (t) - -^-Kl ~ 1)Sjm-2(*), (6a) 

S k i(t) = -S k . ltl+1 (t) - l^S^^t) - ^k(k - l)R k - 2 ,i(t) + ^flV - l)Rk,i-2(t). (6b) 

Note that by definition we have S k o(t) = So k (t) = for all k. We are mainly interested in 
R k i(t), especially Qxif)- It is easy to realize that the equations for the subset {R k i(t)\k+l = K} 
close among themselves provided that the lower-order subset {S k i(t)\k + I = K — 2} has been 
obtained. The equations for the subset {S k i(t)\k + I = K} also close among themselves, 
provided that the lower-order subset {R k i(t)\k + 1 = K — 2} is known. Now for k = I = we 
have Roo(t) = 1 and S 00 (t) = by definition. Similarly for k + l = 1 we have Ria(t) = Roi(t) = 
Sio(t) = Sox(t) = 0. These enable us to solve the cases with k + 1 > 2. 

The case k + I = 2 is of essential interest, because it involves the width of the wave packet. 
By definition we have S 2 o(t) = S 2(t) = 0, and Sn(t) = h/2. These will be useful when the 
case k + I = 4 is considered. For the moment we are interested in the remaining three. It is 
not difficult to find that 



/A; 2 Q 2 (0) + P 2 (0) /iVQ 2 (0)-P 2 (0) n , , Rn(0) ■ . (7 , 

Q 2 (t) = — + — cos2^ + — _ sm2^, (7a) 



^ 2 Q 2 (0) + P 2 (0) ^Q 2 (0) - P 2 (0) . 
P 2 (t) = cos2a;f — /ia;itii(0j sm2o;t, (7b) 

Pn(t) = fln(0) cos2^ - ^ 2 ^°) ~ P2 (°) sin2 ^. (7c) 

2jMjJ 



The width of the wave packet at time t is characterized by the quantity A t x = a/ Q 2 {t), 



and that in the momentum space by A t p = ^JP 2 {t). Thus the width of the wave packet 
oscillates with frequency 2u, just like that for a squeezed ground state. It can be shown that 
both Q 2 (t) and P 2 (t) are positive as they should be. The uncertainty product A t xA t p is also 
oscillating. These conclusions (for a general wave packet) have been discussed in some different 
way previously || ||, so we will not go into more details. A relation that seems not to be 
emphasized in the literature is 

/iV(A t x) 2 + [A t pf = ^u 2 (A x) 2 + (A p) 2 . (8) 



According to this result, A t x reaches its minimum when A t p reaches its maximum and vice 
versa. 
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When k + I = 3, all the Ski(t) can be found to be zero. The equations for 
the Rki{t) can be solved without much difficulty. They are all linear combinations of 
{smut, cosut, sin3wt, cos3cjt}. Since the results are lengthy and not important for further 
discussions we will not write them down. We just point out that they all vanish if the corre- 
sponding initial values Rki(0) are all zero. 

When k + I = 4, it can be found that S 40 (t) = S 04 (t) = 0, S 3 i(t) = 3hQ 2 (t)/2, S 13 (t) = 
3hP2(t)/2, and = 2HRu(t). However, the results for the Rki{t) are rather lengthy, they 

are linear combinations of {1, sm2cut, cos2cut, smAcut, cosAcut}. We only write down one of 
them here: 

_ , . 3^ 4 Q 4 (0) + 3P 4 (0) + 6^Vi? 22 (0) + 3h 2 fi 2 u 2 
W) ~ 8^ 

/i 4 cu 4 g 4 (Q) - p 4 (Q) . , i?i 3 (Q) + ^ 2 i?3i(o) . . 

+ o a a cos2wtH — sm2ut 

fi 4 u 4 Q A (0) + P 4 (0)-Qfi 2 iu 2 R 2 2(0)-3h 2 i2 2 uj 2 A ^ 

+ — — : — : — COs4cut 

8/i 4 CJ 4 

R 13 (0)-fi 2 m 2 R 31 (0) . . . (Q , 
smAut. 9 

We see that higher-order moments about the center have more complicated time dependence, 
but they are essentially oscillating and hence are periodic functions. This is also true for still 
higher ones. Indeed, x t , pt and ip t are all periodic, and so are Rkiif) and Ski(t). 

For larger values of k + l, the solutions are more difficult to find. We will not proceed further 
in this respect. 



III. SIMPLIFIED RESULTS FOR SPECIAL CASES 

The results obtained in Sec. II are rather complicated. In the following sections we will 
discuss the results for special cases with the initial condition 

i> = tp(x, 0) = ip(x - x )e ipox/h , (10a) 
where xq and po are real constants, and (p(x) has definite parity, namely 

ip(-x) = ±<p(x). (10b) 
In this initial state it is not difficult to show that 

x = x , Po=Po, (11) 

and 

i*«(0)=0, 5 fc ,(0)=0, k + l = 1,3,5,... . (12) 

The latter could be derived from the useful relations 

(ipo, (x- x Q ) k (p-p Q ) l i) Q ) = ((fi(x),x k p l (p(x)), 
(^o,(p-PoY(x-x ) k ij ) = (<p(x),p l x k <p(x)), 

A:,/ = 0,1,2,... , (13) 
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which are easy to show. 

The shape of the initial wave packet, \ipo\ 2 , is obviously symmetric about the center xq. This 
symmetry will be kept at later times. To prove this it is sufficient to show that Q2K-i{t) = 
for all natural numbers K. Indeed, we will show that 



R M (t)=0, S M (t) = 0, k + l = 1,3,5,... . (14) 

For k + I = 1 it is true as given in Sec. II. Now suppose that it is true for k + I — 2K — 1, 
and consider the case with k + I = 2K +1. It is easy to realize that the equations for the 
subset {Rki(t)\k + I = 2K + 1} close among themselves and are all homogeneous because all 
of the Ski(t) in the subset {Ski(t)\k + I = 2K — 1} vanish (according to the assumption). 
Since the initial conditions are all homogeneous too, the solutions are obviously Rkiit) = 
{k + l = 2K + 1). Similarly, Ski(t) = for k + l — 2K + 1 as well. The case with k + 1 — 3 has 
been explicitly calculated, and the result is consistent with this general conclusion, as pointed 
out in Sec. II. 

If <p(x) is real, we have another useful consequence: 



R k i(0) = 0, kj = 1,3,5,... . (15) 

As a result the solutions Rki{t) with k + 1 = 2K are also simplified. For example, the sine 
terms in Eqs. (ffa]), ([TED , and (||) all vanish. However, that <p(x) is real is just a sufficient 
condition, not a necessary one for Eq. (|T5|). 

Now if ip(x) in the initial state satisfies the condition 

(<p(x), (xp + px)(p(x)) = 0, ((f(x),/j 2 uj 2 x 2 (p(x)) = ((f(x),p 2 (f(x)), (16) 
we have [cf. Eq. (p~3|)1 

R n (0) = 0, M 2 (0)=? 2 (0), 

and from Eq. (0) we obtain 

Q 2 (f) = Q 2 (0), P 2 (t) = P 2 (0), R n (t) = 0. 

That is, the width of the wave packet keeps unchanged. The second condition in Eq. ( |T6"D 
means that the kinetic energy and the potential energy have the same mean value in (p(x). 
It is easy to realize that all the number states <p n (x) satisfy the conditions ( |10b| ) and (JT6|) . 
However, there exist many other functions that satisfy these conditions. For example, 



<p even (x) = J2ai<p 2ni (x), 0<n!<n 2 <..., ^\ ai \ 2 = 1 (17a) 



i=l i=l 

and 



V odd {x) = J2 b W2n 1+ i{x), 0<m <n 2 < ... , ^|&,| 2 = 1 (17b) 



8=1 i=l 



all have definite parity. And as long as n^+i — rij > 2 for all i, they also satisfy the conditions 
in Eq. (||). 



6 



However, the shapes of such wave packets change with time, though their width keeps 
unchanged. The reason is that Q±(t), and in general the higher-order moments, are still 
oscillating. The conditions for Q±(t), P±(t) etc. to be time independent are 

RiM = i*3i(0) = 0, /iVg 4 (0) = P 4 (0), 2^ 2 Q 4 (0) - 6i? 22 (0) = 3h 2 . (18) 

If tp(x) = (f n (x), a number state, then it can be shown that these conditions are satisfied, as 
expected. However, there exist many other functions that satisfy these conditions. In fact, the 
two classes of functions given in Eq. ( |T7D do if n i+ i — rii>3 for all %. In this case Eq. ( |TBD is of 
course satisfied too. Since now Q±(t) is also time independent, the shapes of the wave packets 
changes with time less than the case where only the width is kept unchanged. In other words, 
their shapes are more rigid. 

In general, the higher-order moments, say Qe{t), of the above wave packets will still be time 
dependent. It will become more and more difficult to discuss the problem in the above manner. 
Since we are not going to find the most general results for these higher-order moments, we will 
proceed in a different way in the next section. 



IV. WAVE PACKETS WITH VARIOUS DEGREES OF RIGIDITY 



In the last section we have found wave packets whose width keeps unchanged with time, 
and ones whose Qxif) up to K = 4 all keep unchanged with time. Obviously the shape of 
the latter changes with time less than the former. Now we introduce a natural number called 
the "degree of rigidity" for a wave packet: If Qxif) is time independent for K — 2, 3, . . . , 2N 
[note that Q (t) = 1 and Q x (t) = 0] but not for K = 2N + 2 (the situation for K = 2N + 1 is 
not important in this definition), we say that the wave packet has degree of rigidity N. Thus 
the two cases mentioned above have degrees of rigidity 1 and 2, respectively. It is obvious that 
the larger N is, the less the shape changes with time. If the width of a wave packet changes 
with time, then its shape changes with time apparently and thus has no rigidity. On the other 
hand, if Qxif) is time independent for all K, then the degree of rigidity for such a wave packet 
is infinity. In other words, its shape is perfectly rigid. Typical examples of such wave packets 
are displaced number states. 

In order to find wave packets with a given degree of rigidity, we consider the quantity 

W kl (t) = (fa, (x - x t ) k (p - p t )%) = R kl (t) + iS kl (t). (19) 
Because ip t = e~ iHt ' Ti "ipQ, we have 



W kl (t) = (^o, (x t - x t ) k ( Pt - p t )Vo), (20) 

where 

Xf = e iHt/h xe- iHt/h = x cos cut + — sin out, (21a) 

IMjJ 



Pt 



iHt/h -iHt/h 



p cos cut — fxcux sin ut. 



(21b) 
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Now we confine ourselves to initial states of the form flTOl). On account of Eqs. (|4]) and (|TT|), 
we have 

Wki(t) = (ipo, [(x — xq) cos tut + — po) sinut] k [(p — p ) cos cut — fiu>(x — x ) sincutjVo)- 

(22) 

Note that for any function f(x) we have (p — po)[e mx ^ h f (x)} = e mx ^ h pf(x), and that the inner 
product is an integral over x, the above equation can be simplified as 

Wki(t) = (tp(x), [xcosujt + (fiu^p sin uit] k (p cos tut — fiivx sin ivt) l (p(x)). (23) 
Using Eqs. ([H]) again this becomes 

W kl (t) = ( Vt (x),x k p l Vt (x)), (24) 

where 

Vt{x)=e- lHt ' h y{x). (25) 

If </?(z) = ^Pn(x), a number state, then (ft(x) = e~ lEnt ' h ip n (x) , and 

W kl (t) = (<p n (x) 1 x k p l <p n (x))=W kl (0), k,l = 0,1,2,... . (26) 

Especially, Qxit) = Wxoit) = Qk{0) for all if, so that the shape of a displaced number state 
is perfectly rigid. This is a well-known conclusion. Here we obtain the conclusion in a different 
and also simple way. 

Consider the function ip even (x) given in Eq. (|17a|) , we have 



¥?r n M = E * exp(-iiW/ftV 2ni (x). (27) 
i=i 

Assuming that — > iV + 1. Because x k p l ip2 ni (x) is a linear combination of 

{(p 2 n i+ k+i(x), (p 2ni +k+i-2(x), ... , <f 2 m-(k+i) (x) } , we have for k + l < 2N that 

CO 

W kl (t) = (^ ven (a;),xVVr n (^)) = H^n^^VW^)) = Wh(0). (28) 

8=1 

Especially 

Qjc(*) = Qjc(0), if = 2, 3,... ,2iV. (29) 

The same result can be obtained for the function ip odd (x) given in Eq. ( |17b| ). [Remember that 
Qxit) = for all odd K as proved in Sec. fH\ ] Thus, if n i+ i — > N + 1, the two classes 
of functions given in Eq. (|T7] ) lead to wave packets [whose initial wave functions are given by 
Eq. ( |10a| )1 that have degree of rigidity not less than N. If some of the differences n i+1 — n,i 
equals N + 1, they will in general have degree of rigidity N. Otherwise the degree of rigidity 
will be greater than N. For large N, the shapes of these wave packets are almost unchanged 
with time. 
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V. SUMMARY 

In this paper we studied the time evolution of a general wave packet in the harmonic 
oscillator potential by examining the time dependence of the various moments Qx{t) of x 
about the center of the wave packet. The differential equations for these objects are derived 
for general K, and are solved for the most general initial conditions up to K = 4. Wave packets 
with constant width, and ones whose Qk{1) up to K = 4 are all constants, are discussed. These 
include displaced number states as special cases. In general the shapes of these wave packets 
still change with time because the higher order moments are oscillating. A natural number 
N, called the degree of rigidity, is introduced to describe qualitatively how much the shape 
is changed with time. The larger N is, the less the shape is changed with time. Displaced 
number states are perfectly rigid and have N = oo. Two classes of wave packets with an 
arbitrarily given degree of rigidity are given explicitly. 
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